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We combine experimental observations by scanning tunneling microscopy (STM) and den¬ 
sity functional theory (DFT) to reveal the most stable edge structures of graphene on Ni(lll) 
as well as the role of stacking-driven activation and suppression of edge reconstruction. De¬ 
pending on the position of the outermost carbon atoms relative to hollow and on-top Ni sites, 
zigzag edges have very different energies. Triangular graphene nanoislands are exclusively 
bound by the more stable zigzag hollow edges. In hexagonal nanoislands, which are con¬ 
strained by geometry to alternate zigzag hollow and zigzag top edges along their perimeter, 
only the hollow edge is stable whereas the top edges spontaneously reconstruct into the (57) 
pentagon-heptagon structure. Atomically-resolved STM images are consistent with either top- 
fcc or top-hep epitaxial stacking of graphene and Ni sites, with the former being favored by 
DFT. Finally, we find that there is a one-to-one relationship between the edge type, graphene 
stacking, and orientation of the graphene islands. 


Introduction 

Edges play a fundamental role in shaping the morphology- - - and electronic- - - properties of 
graphene nanostructures. Electron confinement due to edge boundaries gives rise to energy band 
gaps,—; dI localized states spin-polarization,-^ 12 ' and spin-dependent electron scattering.— Fur¬ 
thermore, graphene edges determine the preferred sites for the attachment of metal atoms,— chem¬ 
ical functionalization, 15 as well as oxygen etching and intercalation. 1617 Understanding and defin¬ 
ing the edge morphology is therefore important to tune the growth of graphene nanostructures as 
well as to modulate the electronic properties of graphene in confined geometries. 

In free-standing graphene, crystallographically oriented edges are of either zigzag (zz) or arm¬ 
chair (ac) type. 6 Graphene with ac and zz edges can be obtained by mechanical exfoliation 1 - and 
etching techniques. 19 However, free-standing unpassivated edges are unstable due to the high den¬ 
sity of dangling bonds, 2 -2r— which induce the zz edges to spontaneously reconstruct into a line 
of pentagon and heptagon pairs, the so called zz(57) or Stone-Wales reconstruction.' 23 -^- A very 
different scenario arises in epitaxially grown graphene, where the interaction with the substrate 
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can stabili z e ac,—zz, 2 -^ 29 and reconstructed edges, 4-30-32 and induce complex graphene-metal 
boundaries 

Despite the potential of epitaxy for tailoring graphene edges evidenced by these studies, pre¬ 
dicting the edge structure is still a challenging task, which requires detailed insight into the inter¬ 
play of substrate interaction and edge morphology. For example, minimal variations in the lattice 
constant and electronic structure of the substrate can lead to very different edge energetics, re¬ 
sulting in hexagonal islands with zz and Klein-type edges in Co(lOOO),- 0 and either hexagonal or 
triangular islands in Ni(l 11), depending on growth parameters. 34-35 The latter strongly differ from 
the quasi-isotropic shape predicted by theory.-^- 

In this work, we show that the experimentally observed morphological transition of graphene 
islands on Ni(l 11) is driven by a substrate-induced reconstruction of zz edges into the zz(57) struc¬ 
ture.This reconstruction, which was found to be unstable on most close-packed metals including 
Ni(l 11) according to previous first principles calculations,- 1 ^, occurs along a particular direction 
with respect to the substrate, and can only be properly reproduced by theory by taking into account 
the stacking of differently oriented edges. 

Methods 

Experiments 

Graphene nanoislands were grown on a Ni(l 11) single-crystal surface kept in a ultra-high vacuum 
(UHV) chamber with a base pressure of 3 x 10 10 mbar. The crystal surface was cleaned by re¬ 
peated cycles of Ar+ sputtering followed by annealing at 800 °C during 1 minute. The surface tem¬ 
perature was measured throughout the experiment using a pyrometer (IMPAC IGA 140). Propene 
was dosed on a freshly prepared Ni(lll) crystal at room temperature. Once the dosing process 
is complete the sample is heated at 500 °C for 5 minutes to nucleate the graphene nanoislands 
and subsequently cooled to room temperature. The islands so obtained have an irregular shape. In 
order to prepare triangular and hexagonal graphene islands (TGIs and HGIs, respectively) we post- 
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anneal the sample to a temperature Ta — 500 °C (650 °C) during a time G — 20 min (10 min). We 
use a propene dose D — 1 L for the TGIs and D — 2 L for the HGIs to compensate for C loss at high 
temperature. The heating rate was 12°C/s and the cooling rate was 1.7°C/s. It should be empha¬ 
sized that the presence of contaminants, such as H or CO, is minimized by annealing the samples 
in UHV. Atomic H as well as molecular species, such as CO, desorb from Ni(l 11) at temperatures 
lower than 300 °C. 36 ' 37 Moreover, we observe that the island shape does not change at 500°C, 
which is another indication that the edge structures are not driven by hydrogenation. Although 
it is notoriously difficult to completely rule out the presence of contaminants, the fact that DFT 
calculations provide a very satisfactory explanation for the STM observations can be considered 
as a further evidence for the absence of contaminants at the edges. More details on the preparation 
of the TGIs and HGIs preparation can be found in Ref.[35j. Topographic images of the surface 
were obtained at room temperature using a variable-temperature scanning tunneling microscope 
(SPECS, STM Aarhus 150) and processed using a freeware software (WSxM 5.0 develop 4.1 ).— 


Theory 

Ab-initio spin-polarized calculations were performed using DFT as implemented in the SIESTA 39 
and ANT.G codes.- 40 For the description of TGIs using SIESTA, we considered a supercell made 
of a three Ni(lll) layers (up to four layers in convergence test, see the Supporting Information) 
and one graphene island containing 22 C atoms, whereas the results using ANT.G were obtained 
employing a graphene island of 33 C atoms placed on a Ni(l 11) surface described by a cluster with 
2 Ni layers. The simulation of HGIs with alternate edge types would require the use of extremely 
large supercells; we have therefore used the nanoribbon geometry to compare the properties of 
different graphene edges on Ni(lll) using SIESTA. The results obtained for the nanoribbon cal¬ 
culations yield information on the edge stability and energetics, which can be extended to the case 
of large hexagonal nanoislands. For such calculations we employed a 4x8 supercell made of three 
Ni(l 11) layers and a graphene ribbon containing 40 C atoms placed on one of the surfaces. In this 
case, only graphene nanoribbons with top-fcc stacking with respect to the underlying Ni(l 11) were 
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Figure 1: Hard sphere models of different zigzag graphene edges on Ni(l 11) for fee (left) and hep 
(right) graphene-Ni stacking. 

considered. Further details on the calculations are given in the Supporting Information. 

Results and discussion 

The small lattice mismatch between graphene and Ni(lll), lxl stacking, and the corresponding 
absence of a Moire pattern make this system an optimal candidate to investigate the structure, 
stability, and epitaxial relationship between the graphene edges and a metallic substrate. Several 
lxl stacking structures have been proposed for graphene monolayers on Ni(l 11) 5 According 
to these previous studies, the top-fee, top-hep and bridge-top configurations are considered to be 
the most stable, with a small energetic preference for the top-fee stacking. Indeed, a recent report 
has confirmed the coexistance of these three graphene configurations on Ni(lll), with a general 
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predominance of top-fee. 46 In the case of graphene nanoislands on Ni(lll) reported here, the 
possible stacking configurations are expected to be the same as for the graphene monolayer. Sym¬ 
metry arguments, however, allow us to exclude bridge-top stacking (see Supporting Information) 
and consider only top-fee and top-hep stacking shown in Figure Q] 

Regarding zz edges, we distinguish two different types according to the position of the outer¬ 
most C atoms with respect to the substrate. We thus define edges having the outer C atoms located 
on top of Ni atoms as zz t and edges having the outer C atoms on hollow sites as zz^, as shown in 
Figure [D The outer C atoms of a zz^ edge can occupy hollow fee or hep sites, depending on the 
stacking configuration and edge orientation. As seen in Figured] the crystallographic directions of 
zzt and zz^ edges with the same stacking differ by 60° and the crystallographic directions of edges 
of the same type but different stacking also differ by 60°. Note also that an edge perpendicular to 
the (121) direction can be either a top-fcc/zz t or a top-hep/zz^ edge, while an edge perpendicular 
to the (112) direction can be either a top-fcc/zzh or a top-hcp/zz t edge. 

In previous work, we have shown that graphene islands that nucleate with random shape upon 
dosing propene on a clean Ni(lll) surface and annealing the substrate above 450 °C evolve to 
either triangular or hexagonal shapes depending on the initial hydrocarbon dose, annealing time, 
and temperature. 35 Such TGIs and HGIs provide an ideal system to study the structure and stacking 
of the edge C atoms on close-packed metal surfaces. Due to the threefold symmetry of the zzt and 
zzh edges introduced above, TGIs present only a single edge type whereas HGIs must alternate 
different edges. By varying the island size and shape it is thus possible not only to characterize the 
different edges but also to study the influence of the edge morphology and graphene stacking on 
the stability of the islands. 

Figure (2fa) shows a representative STM image of TGIs obtained by annealing the substrate 
to 500 °C and subsequent cooling to room temperature (see Methods). Under these conditions, 
most islands exhibit straight edges, which indicates that they can only be of either zz or ac type. 
Topographic images with atomic resolution on both graphene and Ni are shown in Figure |2tb) 
and (c). We observe that the TGIs have unreconstructed zz edges with no defects, in agreement 
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with previous results. 34,35 The edge shown in Figure [2[c) runs perpendicular to the [112] direction, 
i.e., parallel to a high symmetry [110] direction of the Ni lattice. A honeycomb and an hexago¬ 
nal lattice, representing the atomic positions of the graphene and surface Ni atoms respectively, 
have been superposed on the image for illustrative purposes. Note that the STM image taken over 
the graphene-covered region exhibits almost hexagonal contrast, despite the honeycomb atomic 
structure of graphene. This is a well-known electronic effect induced by the graphene-metal in¬ 
teraction, which breaks the symmetry between C atoms in on-top and hollow positions. 46-48 In 
particular, previous theoretical simulations concluded that in the case of graphene on Ni(lll) the 
bright spots in the STM images at zero or very small bias correspond to C atoms located in hollow 
positions,^ 47 as reproduced also by our calculations. Besides, all the images shown in this work 
have been acquired under well controlled tip conditions, which allows us to exclude tip-induced 
contrast effects. All this suggests that, away from local variations related, e.g., to surface impu¬ 
rities, lattice imperfections, and the proximity of the edge, the spots with larger intensity in the 
STM images correlate with C atoms in hollow positions. In any case, the hexagonal units of the 
honeycomb lattice can be clearly identified and extrapolated to the Ni lattice, which let us draw 
two important conclusions: i) that vertices of the honeycomb lattice at the graphene-Ni lateral 
boundary (red circles) correspond to hollow positions on the hexagonal lattice, and hence the edge 
is of zzh type; ii) The vertices representing the other sublattice (blue circles) are located on top 
of the Ni atoms, and hence graphene has either top-fee or top-hep stacking with the substrate, in 
agreement with our initial hypothesis. 

The great majority of the TGIs in our samples point in the same direction, as seen in Figure[2ta). 
Because there is a one-to-one relationship between the edge type, stacking, and orientation of TGIs 
(Figure S2 of the Supporting Information), this implies that the TGIs have a preference for a unique 
edge type and stacking combination, which, according to our analysis, corresponds to zz^ and 
either top-fee or top-hep. Moreover, we find evidence that the edge type dominates the preference 
for the epitaxial stacking of the inner C atoms in the islands. The enlarged area of Figure |2[a) 
shows two TGIs pointing in opposite directions, meaning that they are rotated with respect to each 
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other by 60°. This inverted orientation can only be explained by a change in the stacking or in 
the edge type (see Figure S2 of the Supporting Information). Note that all the TGIs pointing in 
the direction opposite to the predominant one [indicated by arrows in Figure Eta)] have a small 
size, below 10 nm 2 . Since the effect of the edge type on the system’s energy is higher for smaller 
islands, the most plausible scenario is that the TGIs with inverted orientation maintain the lowest 
energy edge configuration and change their stacking. 

These experimental observations have been confirmed and complemented by DFT calculations. 
Figure [3] shows the model systems used in the simulations, namely TGIs with zzt and zz^ edges 
and top-fee or top-hep stacking. The TGI and the top two Ni layers of a three layer Ni slab are 
fully relaxed, which leads to significant plastic deformation of the islands owing to edge-induced 
relaxations. These effects are expected to be significant in small size islands such as the one used 
here, where the size is imposed by computational limitations. However, results obtained on larger 
islands using a slab that contains only two metal planes confirm the edge energies obtained for the 
smaller TGIs (see Supporting Information). The relative energies of different edge types are cal¬ 
culated by subtracting the total energy of the TGIs shown in Figure Eland dividing by the number 
of C atoms. This method neglects the plastic deformation of the islands, which would be difficult 
to disentangle by itself and that we estimate to have smaller weight compared to edge effects. Our 
results indicate an energy difference of ~ 0.67 eV between C-edge atoms in hollow and top posi¬ 
tions. From this, we estimate an energy difference between zz^ and zzt edges of A E zz ~ 0.27eV/A 
(see the Supporting Information for more details). This value might be overestimated because 
different edges induce different relaxations in the inner part of the islands. However, this effect 
should be mitigated for the calculations of larger islands reported in the Supporting Information. 
In the following, we use the average value of all these calculations to set the difference of the zz^ 
and zzt edge energies, while our maximum and minimum values are used to associate an error bar 
with this result. According to this estimation A E zz = 0.22 ±0.05 eY/A. Therefore, DFT predicts 
that the most stable TGIs have zz^ edges, in agreement with the experimental results. The prefer- 
ence for zz^ edges can be understood by considering the deformation of the sp hybridized orbitals 
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of the edge C atoms as they bind to the surface Ni atoms. In the case of zz^ edges only a small 
deformation is required to form a bond with a Ni atom, whereas for zz t edges a larger deformation 
takes place, which also leads to the bent edge structure that can be observed in Figure [3] 

With respect to the stacking with the substrate, the calculations predict a small energy gain for 
the top-fee stacking with respect to the top-hep, of the order of a few tens of meV per C atom, 
which is not very significant in the case of the small islands considered in this study. This result 
agrees with the calculations for extended graphene on Ni(lll). 4j Moreover, it is consistent with 
the preferential orientation found for the TGIs. The small energy difference obtained for the two 
stacking configurations is reflected by the fact that islands with opposite orientation, although very 
rare, can still be observed when their size is small enough [Figure [2j a)]. 

Figurea) shows that TGIs evolve into HGIs by annealing the Ni(lll) surface up to 650 °C 
for 10 minutes. This seems to imply that the growth of the TGIs is kinetically limited, whereas 
the equilibrium shape of the graphene islands on Ni(lll) is hexagonal, as confirmed by the re¬ 
sults of our DFT simulations. This is similar to the homo- and hetero-epitaxy of metals on (111) 
surfaces.- 49,50 The triangular and hexagonal island shapes allow us to compare the structure and 
stability of different edge types. Contrary to TGIs, HGIs are constrained to exhibit zz^ and zzt 
edges alternated along their perimeter. Because of the 60° angle between adjacent edges, the for¬ 
mation of ac edges can be excluded since these would be oriented at 30° with respect to the zz 
edges. As in the case of TGIs, HGIs have their edges aligned with the high symmetry directions of 
the substrate, and hence possess only zz edges. Since zz^ and zz t edges alternate independently of 
the stacking (Figure H]), the edge energy contribution is equivalent for top-fee and top-hep stack¬ 
ing. The orientation of the islands does not change as they evolve from triangular to hexagonal 
by increasing the annealing temperature from 500 °C to 650 °C, which indicates that the prefer¬ 
ential stacking for the HGIs and TGIs is the same, that is, top-fcc according to the DFT results. 
Figure |4[b) shows a detail of an HGI. Edges 1, 3, and 5 appear atomically straight and correspond 
to zzh edges. Edges 2, 4 and 6, on the other hand, present a few structural imperfections and are 
shorter compared to the odd-numbered edges. These edges should be of the zz t type according 
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to the hard-sphere models presented in Figure Q] However, atomic resolution images of the HGI 
show that these edges undergo a reconstruction that has a periodicity of two benzene rings. This is 
clearly visible in Figure BJc), where an island corner is shown at the intersection of edges 1 and 2. 
We assign this reconstruction to a pentagon-heptagon zz(57) structure, which has a periodicity of 
two benzene rings and matches the observed atomic structure. 

In order to confirm these results, we have performed extensive DFT calculations of graphene 
nanoribbons on Ni(lll). Ribbons allow us to study two edge types simultaneoulsy, as HGIs, but 
are computationally less demanding compared to hexagonal islands. In the following we only con¬ 
sider the zz(57) reconstruction, since the other known 2x2 reconstruction of a zz edge, labeled 
as zz(ad) in Ref. |3li is energetically excluded by calculations shown in the Supporting Informa¬ 
tion. Figure [5] presents the optimized structure of graphene nanoribbons with zz^, zzj 1 (57), zz t and 
zz t (57) edges on Ni(l 11). For such structures, as described in detail in the Supporting Information, 
we obtain the following relationships between the energy of the four edge types: 


E zzt(5i) — E zzt 0.15eV/A 

E zzh( 57) = E zzh + 0.16 eV/A (1) 

E zzh(51) + E zzt(57) = 1-07 eV/A. 


These numbers show that it is energetically favourable for zz t edges to undergo the 57-reconstruction 
( E zzt{5l) < E zzt), whereas the opposite is true for zz^ edges (E :zh{57} > E zz h), in agreement with the 
experimental observations. 

From our calculations of the energetics of adsorbed nanoribbons it is only possible to obtain, 
once stacking is taken into account, three relations (Eq. [U for four unknowns. Therefore, addi¬ 
tional information is required to estimate the energy of the different edge types. Here we use our 
estimation of A E zz obtained from the calculations of the TGIs. In this way we obtain the edge 
energies represented in Figure [bja), where the edges with lowest energy are indeed the zz^ and 
zz t (57) types. 
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The results of our calculations are in disagreement with recent theoretical studies of graphene 
nanoribbons on Ni(l 11),-— in which the unreconstructed edge is predicted to be more stable than 
the reconstructed one. However, these studies do not consider the influence of stacking, which 
leads to the conclusion that the 57-reconstructed edges are always less stable than the unrecon¬ 
structed ones. Besides, at least in one of these works/ the estimation of the edge formation 
energies is based on the assumption that the ribbons are perfectly symmetric, which is not true 
if epitaxial stacking is taken into account (see Figure S4 of the Supporting Information). However, 
it is interesting to note that the average value (over the two stacking configurations) of our edge 
energies is in good agreement with the values reported by Gao et al . 31 using plane-wave calcula¬ 
tions and the same exchange-correlation functional that we use here (Perdew, Burke and Ernzerhof 
functional 5 —, see the Supporting Information). The agreement is particularly good for the unrecon¬ 
structed zigzag edges, for which our average edge energy is 0.53 eV/A, exactly the value reported 
by Gao et al. for the zigzag edge energy (without resolving the stacking dependency). Our average 
value for the 57-reconstructed edge energies is ~0.54 eV/A, which is somewhat lower than the 
0.60 eV/A reported by Gao et al., but still in good correspondence. 

With the edge energies in Figure [Ha), following the classical Wulff construction, the equilib¬ 
rium shape of graphene islands on Ni(l 11) turns out to be hexagonal, as shown in Figure (Hb). The 
different length of the zz^ and zz t (57) island edges observed in Figure 0Jb) can also be accounted 
for by such a model, and reflect the larger stability of zz/ z edges as compared to zz f (57) edges. On 
the other hand, if no reconstruction would occur for the zz t edge, the islands would have an almost 
triangular shape, as shown in Figure (Hb) (the optimal shape calculated for triangular islands ex¬ 
hibits also portions of reconstructed armchair instead of zz t edges, see the Supporting Information 
for more details). 

These results outline a possible scenario to explain the evolution of TGIs into HGIs on Ni(l 11) 
when the annealing temperature is increased from ~500 °C to ~650 °C. According to our theo¬ 
retical model, HGIs with alternate zz^ and zz,(57) edges have lower energy compared to TGIs of 
equal size (we estimate HGIs of 10 nm 2 to be ~3 eV more stable than TGIs of the same size, see 
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the Supporting Information). This suggests that an energy barrier must be overcome in order to 
achieve the equilibrium hexagonal shape. In view of our results, we can speculate that this energy 
barrier is associated with the 57-reconstruction of the zz edges. Below a certain temperature the 
system does not have enough energy to overcome this barrier and, as a result, the zzt edge does not 
reconstruct and the island grows into a triangle. On the contrary, as the temperature increases the 
barrier can be easily surpassed and the zz t edges efficiently reconstruct into zz t (57), giving rise to 
hexagonal nanoislands. 

Conclusions 

In summary, we have shown that the structure and stability of the edges of graphene islands grown 
on Ni(l 11) are dominated by the stacking of the edge atoms relative to the substrate. Conversely, 
the edge type and energetics determine the island shape and, in TGIs smaller than 10 nm 2 , also the 
stacking relationship of the inner C atoms with the substrate. Atomically-resolved STM images 
show that TGIs and HGIs exhibit only zz-like edges. As the epitaxial constraint imposed by the 
Ni substrate breaks the six-fold symmetry of free-standing graphene, we distinguish between zz^ 
and zzt edges, which differ in the position of the outermost C atoms relative to the substrate lattice. 
TGIs are bound uniquely by zz^ edges, whereas HGIs are bound by alternate zz^ and reconstructed 
zz t (57) edges. Accordingly, DFT calculations show that the energy of the zz^ and zz t (57) edges 
is about 0.2 eV/A smaller relative to the unreconstructed zzt edges, which are not stable on this 
surface. The edge energetics fully accounts for the shape of the HGIs observed experimentally 
and suggests that the temperature driven transition from TGIs to HGIs is an activated process 
related to the existence of an energy barrier for the 57 pentagon-heptagon reconstruction. The 
stacking of the TGIs and HGIs is experimentally determined to be either top-fee or top-hep. All 
the TGIs and HGIs larger than a few nm 2 adopt the same stacking, which, according to DFT, 
corresponds to the top-fee configuration. The above considerations are important when a well- 
defined lxl stacking leads to strongly orientation-dependent edge stacking configuration, as is 
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the case of Ni(l 11) or Co(0001).- 0 This can apply in a different degree also to other close-packed 
metal surfaces. Likewise, the stability of zz edges and their tendency to reconstruct may play a role 
in determining the density and orientation of grain boundaries in extended graphene layers grown 
on metal substrates. 
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Figure 2: (a) STM image of a large sample area with TGIs of different sizes. The enlarged image 
shows two small TGIs oriented in opposite directions. Setpoint current and bias voltage: I — 1.9 
nA, Vt — 3.7 mV. (b) STM image of a single TGI. Setpoint: 7=1.0 nA, Vb — 0.9 mV. (c) Zoom- 
in of an edge region with the image contrast set to show the atomic resolution on both graphene 
and Ni layers. The left part of the image corresponds to the graphene edge with a superposed 
honeycomb lattice (red and blue circles), the ri^fit part to the Ni(lll) surface with a superposed 
hexagonal lattice (green mesh). The black stripe is a transition region related to the finite size of 
the tip where atomic resolution is missing. 
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Figure 3: Relaxed model configurations of TGIs on Ni(l 11) calculated by DFT, together with the 
corresponding energies per C atom in meV. The latter are given with respect to Eo, the energy of 
the most stable TGI. 
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Figure 4: (a) Large area STM image of HGIs. Setpoint: I — 1.9 nA, Vb — 2.1 mV. (b) STM image 
of a single HGI. Setpoint: / = 9.6 nA, V/, = 0.9 mV. (c) Detail of a comer of the HGI shown in (b) 
with zzh and zzt(57) edges. The black lines in (c) illustrate the doubling of the periodicity observed 
in the reconstructed edge. 
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Figure 5: Relaxed edge structures calculated for graphene nanoribbons with a width of 4 benzene 
rings on Ni(l 11). For visualization purposes, the edge C atoms are represented in green. 
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Figure 6: (a) Formation energies of graphene zz edges on Ni(lll). (b) Equilibrium shape of 
graphene nanoislands on Ni(l 11) obtained by minimizing the edge energy. If edge reconstruction 
is allowed the islands exhibit zz/, and reconstructed zz f (57) edges, adopting an hexagon-like shape. 
If reconstruction is inhibited the islands have a triangle-like shape and exhibit zz/, and reconstructed 
armchair edges. 
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